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. Abstract 

O 

Solutions of the Polchinski exact renormalization group equation in the scalar 
O(N) theory are studied. Families of regular solutions are found and their rela- 
^h. tion with fixed points of the theory is established. Special attention is devoted 

to the limit N = oo, where many properties can be analyzed analytically. 
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1 Introduction 

The Exact Renormalization Group (ERG) has proven to be a powerful non-perturba- 
tive method for studies of phenomena in quantum field theories at various scales. 
It stems from the Wilson renormalization group [|] (see, for example, Ref. || for a 
review) and is, basically, its continuous version adapted to quantum field theory ||-|| 
(see Refs. 0-|[n|] for reviews). The ERG method allows for general non-perturbative 
studies of low-energy effective actions, in particular their flows and fixed points, and 
is also recognized as a rather powerful calculational framework. 

The central object within the ERG approach is the (Wilsonian) effective action 
S*a or the Legendre effective action, which gives an adequate description of the theory 
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at the scale A in terms of renormalized quantities. The action satisfies an equa- 
tion in functional derivatives which determines completely the renormalization group 
evolution of S\ with the scale. This equation is called the ERG equation. 

Various formulations of the ERG have been proposed, the most widely used are 
the Wegner- Houghton sharp cutoff equation ||, Polchinski ERG equation for the 
Wilson effective action (see Ref. [12j for detailed derivation and discussion), and 



equations for the Legendre effective action (|T^]-[[T7|. All these equations are essentially 
equivalent, the relation between them was clarified in Refs. |Tj]] and ||18||. It turns out 



that in practice they are too complicated to be solved exactly. Two well-developed 
non-perturbative approximate schemes for solving the ERG equations are: (1) trun- 
cation in the operator content to a finite set of operators and (2) expansion in 
powers of space-time derivatives p0 |- ||22|| . Whereas the first scheme is of limited relia- 
bility 



23, 24 



the second one, called the derivative expansion, was shown to be quite 
reliable and accurate enough in finding the fixed points and calculating the critical 
exponents || [H], |22|,[p^|-|07|. In particular, the fixed point (FP) corresponding to 



the Ising model in three dimensions || ftB|],P5|-[37|, and the conformal FPs in two 
dimensions p7 , [32] were calculated using the derivative expansion. The lowest order 
of the derivative expansion, O (d°) -order, also called the local potential approximation 
(LPA), already gives rather good accuracy in a large variety of problems (see Refs. 
and ||25|| ). The next-to-leading order (0(<9 2 )-order) usually improves the results of the 
leading order and has also been thoroughly studied in the literature (see, for example, 
Refs. m, H, m and |HJ). 

To outline some general features of the ERG approach let us consider the case of a 
scalar field in d dimensions. A generalization to the fermionic case is straightforward 

m 



the gauge invariant formalism is, however, more involved ||39|| . For a scalar field 
(f) the effective interactions within the LPA are restricted to a potential Va(</>) of a 
general form and do not include the derivatives. The ERG flow equation reduces 
to a non-linear second order partial differential equation which symbolically can be 
written as 

/(0, t) = T{f"(<j>, t), /'(0, t), /(</>, f); v, d), (1) 
where t = — ln(A/Ao) is the flow parameter, Aq is a scale where initial conditions are 



set, 



dV A {4>) 



and r] is the anomalous dimension. The prime denotes derivation with respect to the 
field 4> and the dot denotes derivation with respect to the parameter t. The concrete 
form of the ERG equation will be given in Sect. g. 

FP solutions /*(</>) satisfy Eq. ([[]) with / = 0, i.e. the equation 



W'(0),/'(0),/(0);r ? ,rf) = O. 



(2) 



For the LPA to be consistent the value of the anomalous dimension in the FP equation 
must be set to ^ = 0. Eq. (fj) is supplied with two boundary conditions, usually at 
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= 0. Assuming that the theory is symmetric under ^-transformations <fi — > —<p, 
these conditions are often chosen as 

/(0) = 0, f(0)= 7 . (3) 

A consistent ERG FP equation always has the trivial solution /* = 0. It describes the 
free massless theory and is called the Gaussian FP (GFP). For all known formulations 
the ERG equation is stiff and most of its solutions singular at some finite value of 
the field, making them not acceptable from the physical point of view. Moreover, it 
turns out that only finitely many solutions of Eq. (0) do not end up in a singularity, 
thus giving massless continuum limits with the prescribed field content f22fl . They 

correspond to particular values 7 = 7* , n = 1, 2, 

Let us consider Eq. (0) in a wider mathematical context and study more general 
solutions, namely solutions for arbitrary 77. They are parameterized by 7 and 77 
(for d fixed), i.e. / = /(</>; 7, 77). The dependence of T on 77 is smooth enough, 
as a consequence /(0;7, 77) depends at least continuously on 7 and 77. It is easy to 
understand that regular solutions correspond to certain points (7,77) in the space of 
parameters which form continuous curves 77 = 77^(7). The special values 'yi^, which 
were already introduced above, satisfy 7/ n (7* n ^) = 0. In what follows, by abuse of 
terminology, we will be referring to these regular solutions of Eqs. (§) and (|3|) with 
arbitrary 77 as FP solutions, keeping in mind that physical FPs in the LPA correspond 



to 77 = 0. Such curves of regular solutions were first found in Ref. |32 . 

Obviously, it is useful to know the complete space of solutions of the ERG equation 
in the leading approximation. There is also a practical reason for studying FP problem 
(@); (0) with 77 7^ 0. The FP equations to next-to-leading order, i.e. the equations 
obtained with both O(d ) and 0(d 2 ) terms in S taken into account, form a system 
of coupled nonlinear differential equations. They are stiff and solving them require a 
simultaneous fine tuning of 77 and 7. This makes the direct integration of the system 
a hard task, and one of the ways to solve it is to use an iterative numerical procedure 
2B|. In some cases, as for example for d — 3, one can perform such iterations starting 



with a solution of the first order FP equation with 77 = 0, i.e. with a physical FP.[] 
However, for d = 2 this is no longer possible. It turns out that due to the nature 
of the problem only periodic or singular solutions, neither of which is a physical FP, 
exist for 77 = |]27|, ^ (see Sect. 2 for details). For this reason one needs leading 
order FP solutions with rj 7^ as initial solutions to start the iteration procedure. 

In the present article we study solutions of the Polchinski ERG equation in O(N) 
scalar field theories. Though much research on the structure of the space of solutions 
has been already done, we feel that some questions remain unanswered. Namely, in 
all previous articles the condition of regularity was imposed in the course of numerical 



integration of the equation (see Refs. p^, 23. E6| and [29] for algorithms). No "an- 



alytical" understanding of how this condition works and how it gives rise to critical 



This solution is the Wilson-Fischer FP for d = 3. 
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curves 77^(7) was gained. Also, the relation between the curves and the asymptotics of 
the solutions /(</>) for large remains obscure. To explain this last issue we note that 
solutions of the Polchinski ERG equation in the LPA have, for large 0, the asymptotic 



form 26 



M 7, V) ~ (l - I) <t> + D(v, + ■■■, (4) 



where we introduced 



d 77 



The coefficient .0(77,7) in expansion (f|) cannot be derived from the asymptotic anal- 
ysis of Eq. (H) alone, without solving the FP equation. Fixing D is equivalent to 
setting the boundary conditions at <fi = 00, and only a countable number of values of 
this coefficient correspond to FP solutions. The values depend on 7 and r] n in some 
intrinsic way. In the present article we will address these and other related issues. 

We found it convenient to study them in the context of O(N) scalar theories. The 
reason is that for N = 00 a general solution of the lowest order ERG flow equation 
is known analytically though indirectly [^, [5| |41| . Namely, the inverse function 
= 4>{f, t) can be obtained in a closed analytical form. This will be enough for 
our purpose. We will study various properties of the corresponding FP solutions for 
N = 00 and relate them to the analogous properties of FP solutions for finite N, 
where only numerical or approximate analytical results are available. In particular, 
we will analyze the appearance of the critical curves 77 n (7) of regular solutions for 
N = 00 and see how they change for finite N. 

We would like to note that O(N) scalar models (also called spherical models) are 
interesting on their own |4l|. They play an important role in the description of phase 
transitions. In particular, N = 4 is relevant for the QCD phase transition |2|, |14|, [EJ. 
These models were also studied within the ERG using the 1/iV-expansion [Q (see 
Refs. |l(J and |30| for other studies and references therein). 



The plan of the article is the following. In Sect, g we introduce basic notations and 
discuss the Polchinski ERG equation in the O(N) scalar model for finite N. Sect. [3] 
is devoted to the Polchinski ERG equation in the limit N = 00. We formulate the 
condition of regularity of solutions, find regular solutions and study their properties 
in detail. The appearance of critical curves rj = 77 n (7) is analyzed. In Sect. || we 
analyze the behavior of the critical curves for finite iV and show that they match the 
N = 00 case. Sect. |5] contains some discussion and concluding remarks. 



2 0{N) Polchinski ERG equation 

In this article we consider a general 0(iV)-symmetric scalar field theory in <i-dimen- 
sional Euclidean space. 
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Let us denote the field components by <p a (x), a = 1, 2, . . . , N and assume that the 
Wilson effective action at the scale A is of the following general form: 

SaM = \J " M(p 2 , A) • 0_ p + S int M- (6) 

Here 0^ denotes the field in the momentum representation, i.e. the Fourier transform 
of (j> a (x). In fact the field depends on the scale A as well, i.e. 0^ = (p a (p; A). The dots 
in Eq. (|6|) stand for the contraction over indices of O(N). The inverse propagator is 

M ab (p 2 ,A) = 5 ab 1 



iV,A) : 



K(p 2 /A 

p2 



2\ 



where K{p 2 /A 2 ) is a cutoff function (also called a regulating function or simply a 
regulator) || |TIJ. It has to satisfy the following two requirements: (1) K(0) = 1, 
and (2) K(z) — > as z — > oo fast enough. The second property guarantees that 
contributions of high momentum modes are suppressed. In some studies on the ERG 
a sharp cutoff function, defined by K(z) = 1 for < z < 1 and K{z) = for z > 1, 
was used || |6|. This choice corresponds to integrating out modes with p 2 > A 2 
exactly. In the present article we assume that K(z) is a smooth enough function, so 
that higher momentum modes are integrated out effectively. The term Si n t t A[4>] is the 
interaction term. In general it contains all possible 0(AQ-invariant interactions. 



2.1 Polchinski ERG equation for arbitrary TV 

The central idea of the ERG approach is essentially the following: under a change of 
the scale A the Wilson effective action changes in such a way so that the generating 
functional of the theory and, therefore, the Green functions remain unchanged || |J. 
This requirement gives rise to an equation on S\[<f>] which is called an ERG flow 
equation. One of the most widely used versions is the Polchinski equation: 




5S 



5ip- p 5ip p 



(P ,A)^ P • 



5S 

5ip p 



(7) 



where we introduced = 0*/vJV and S = S[(p;A] = (1/N)Sa[(/)], the notations 
which will be convenient in the N — ► oo limit. This equation for N = 1 was first 
derived and discussed in Ref. H, for TV > 1 it was studied in Ref. 



29 



The structure of Eq. (|7]) becomes more transparent if it is rewritten in terms of the 
dimensionless momentum p = pj A and the dimensionless field variable = A 1+d / 2 ^. 
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We also introduce the renormalization group flow parameter t = — ln(A/Ao), where 
A is some fixed ultraviolet reference scale. When expanded in powers of fit a general 
action of interaction has the form 



oo „ 

Smt[<P',t] = J 



k=l 
ai 

Pi ~ r P2 



d d pi . . . d d p 2k 

(27T) d ( 2 



x <p%<p°* 



k-i) 



■n;ai,...a 2 k (Pi, • • -,P2k]t) 



Pl + ... +p2k>- 



Note that the terms of the expansion contain only an even number of fields and are 
0(A^)-invariant. The Polchinski ERG equation, Eq. (|7|), becomes 



as 

~dt 



d d p 



(27r) d 



+ Sd + 



k'(p 2 : 

d d p 



5S 5S 



S0-p 

d 

1 

2 



80p 

T,{t) 



;Tr 



6 2 S 



N \50-p50p t 
,K'(f) 



{2n) d 
d d p A „ d' 5S 



K(p 2 



5S 
80 p 



(9) 



The partial derivative in the l.h.s. acts on the t-dependent coefficient functions 
s n ;ai,...a 2 h m ^q. @ and does not act on 0*. The scale dependence of the field 0® is 



determined by 



d0 



dt 



d 7] 
- + - 
2 2 



<P 



where r] is the anomalous dimension, and has already been taken into account in Eq. 
(P) by adding the rescaling terms (the second line in Eq. (|9])). The prime in the 
momentum derivative in the last term in the r.h.s. means that it does not act on the 
5-f unction of the energy- momentum conservation (see Eq. (Q)). 

Being supplied with an initial condition S[0; t = 0} = S o [0] at t = (or A = A ), 
where So [0] is some given functional, Eq. @ determines the flow of the effective action 
and allows to calculate, at least in principle, the action and physical characteristics 
of the theory at low energies, when A — > (or t — ► +oo). The limit Ao — > oo, when 
the cutoff scale Ao is removed, is called the continuum limit and corresponds to the 
renormalized theory. 



2.2 Leading order approximation 

As was mentioned already in the Introduction, to find solutions of the ERG equation 
in concrete calculations various approximate (but non-perturbative) techniques are 
used. A reliable and efficient one is the derivative expansion. In the coordinate 
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representation it corresponds to the following expansion of the action of interaction 
in powers of derivatives: 

Sint[<P',t] = 
+ 

where the potentials V, Z and W do not contain derivatives of the field. We also took 
into account that these functions depend on the field only through the 0(iV)-invariant 
combination (p 2 = (p[x) ■ 0(x) = a (x)ip a (x). Restricting the action to the first term 
gives the leading-order of the derivative expansion, the LPA. In this approximation 
the Polchinski ERG equation reduces to a partial differential equation for the function 
V(z 2 ;t), where we denoted (p 2 = z 2 . Introducing the notation f(z,t) = ■^V{z 2 ]t) 1 
one can write the equation as 




/ = If - 2ff + 



1-1U + A-* 

NJ z 



f' + 



A H 



1 



1 ) 1 

NJ z 2 



f, 



(10) 



where the dot and prime denote the derivative with respect to the flow parameter t 
and variable z respectively and, as before, 



2 2 



(see Eq. (|5])). FP solutions f(z) are solutions of Eq. (p~0|) with the property / = 0, 
i.e. they satisfy the equation 



N 



2ff + 



1 



1 \ 1 . 
— - + A" 

N z 



f' + 



A^ 



i4U 



NJ z 2 



f = 0- 



'IV 



Note that rj is constant here. 

In this paper we will assume that the potential V(z 2 ) is analytic in z 2 at z 2 = 0, i.e. 
V(0) is finite and its expansion at z = is convergent and contains only non-negative 
integer powers of z 2 = <f) 2 . This assumption is fulfilled in most of scalar models and, 
therefore, does not add new essential restrictions in field theory applications of the 
ERG. Correspondingly, the expansion of f(z) at z = includes only odd powers of 
z. To have a unique FP solution we must add two conditions on f(z). In field theory 
problems it is natural to impose them at z = in the form 



f(z)\ z=o = 0, f'(z)\ z=0 = r (12) 

For purely historical reasons these conditions are called initial conditions. The first 
condition is a consequence of the 0(iV)-symmetry of the potential V(0 2 ). The pa- 
rameter 7 plays the role of the square of the mass of the field. 
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For further studies it is convenient to define a new variable y and a new function 
u(y,t) instead of z and f(z,t) in the following way: 

z = y/ty, (13) 

/CM) = ^f«(y,f). (14) 

Note that these formulas define u(y,t) on the semi- axes y > only. In terms of y and 
u(y,t) Eq. ( |TD[ ) becomes 

«(y, *) = |<„ - %K - M ' + (i + § + 2 f A ~) K + su > (1 5 ) 

where we introduced the notation 

s = A + + A~ = 2-r). (16) 

The prime in Eq. (|j|) denotes the derivative with respect to y. The FP equation 
now reads 

- 2yuu' y -u 2 +(l + ^ + 2yA") < + SM = 0. (17) 

Let us rewrite initial conditions (|T2D in terms of u(y). Using relations flT3|), (|14D it 
is easy to check that the condition /(0) = together with the assumption about the 
behavior of f(z) at z = are equivalent to the requirement that u(y) is analytic at 
y = 0, i.e. w(0) is finite and its expansion at y — is convergent and contains only 
non-negative integer powers of y. Therefore, from now on we will be interested only 
in solutions of Eq. (|17]) which are analytic at y — 0. The second initial condition in 
( |T2| ) transforms to 

u(y)\ y=0 = 2 1 . (18) 

Using the analyticity of u(y) the value of u'(y) at the origin can be calculated from 
the FP equation 

2 7 (2 7 - s) 
u {y)\y=o = , 2_ ■ 

N 

In fact, to be physical solutions functions f(z) or u(y) must satisfy even stronger 
requirements. Recall that f(z) is the derivative of the effective potential V with 
respect to the field variable. In field theory applications one is interested in solutions 
such that V(0 2 ) and its derivatives are finite and do not have singularities for all 
finite values of the field variable. This means that f(z) does not have singularities 
for finite z or, equivalently, u(y) does not have singularities for finite y > 0. Such 
solutions will be called regular solutions. 
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From Eqs. ([TT]) and fll2"|) , or, equivalently, from Eqs. (|T7|) and ([18]), one can see 
that for a given c? a FP solution is completely determined by the values of 7 and 77. 
Correspondingly, any solution of the FP problem can be represented by a point (7, 77) 
in the parameter space. 

One comment is in order. Since the potentials Z and W are not taken into account 
in the LPA, the kinetic term in the effective action, which contains two derivatives 
of the field, matches ERG equation (|S|) for r\ = only. Regular solutions of Eqs. 
(|TT|), (|TT|) with 77 = will be called physical FPs. As we have already explained in 
the Introduction, we will relax this condition and study a wider class of FP solutions 
with 77 7^ 0. To summarize, in the present article we will study solutions u(y) of Eq. 
(|T7| ) with arbitrary 77, both regular and non-regular, which satisfy condition ([18]) and 
are analytic at y = 0. 

The asymptotics of a solution u(y) for large y can be obtained from ERG equation 
([IT]) by a straightforward calculation. Namely, one writes a general expression for the 
expansion of u(y) as a series in decreasing powers of y and fixes the coefficients and 
the exponents from Eq. flT7|). One gets 

, , n aD 2 r, R D 3 a(a + 1) n R l + %3 a , . . 

u(y) = s-Dy 13 - — -y w - — / > y w - D— + ■ ■ • , 19 

vy; y s(2-a) y s 2 (2 -a) 2 s(2-a) 2 ' V ; 

where 

_ 2A" 2A~ 1 s 

a ~ ~s~ ~ A+ + A"' ~2-a ~ _ 2A+' 

As was already mentioned in the Introduction, this technique does not allow to fix 
the constant D = .0(7,77). Expansion ([19]) is, of course, in a full correspondence with 
asymptotic formula @) with 




The FP equation has two trivial solutions. The first one is obtained for 7 = and 
is equal to / = 0. This is the GFP describing the free massless theory. It corresponds 
to the line 7 = on the (7, 77)-plane. The second solution is obtained for 7 = s/2 
and is given by the linear function f(z) = sz/2, where s is defined in Eq. ([161) . These 
solutions, termed high temperature FPs or infinite mass GFPs, are described by the 
straight line rj = 2 — 27 on the (7, r/)-plane. For the sake of brevity we will call them 
trivial FPs (TFP). As one can easily see from Eqs. ( |Tj] ) and (0), these solutions 
correspond to u = and u = s respectively. 

2.3 N = 1 case 

Before starting the analysis of FP solutions in the O(N) model for arbitrary N it is 
instructive to consider the N = 1 model. In this case the 0(A^)-symmetry reduces to 
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the ^-symmetry under the reflection <p — ► — 0. The leading order Polchinski ERG 
equation (|TT|) with initial condition (|l^) becomes the familiar FP problem (see |26] ) 



f"(z)-2f(z)f'(z) + A-zf(z)+A + f(z) = 0, (20) 
/(0) = 0, /'(0)= 7 - (21) 

For A~ = or A + = this equation can be integrated analytically. Let us analyze 
the case A - = first. With this condition, which implies A + = d, Eqs. (20) and ([Hp 
become 

f(z)-2f(z)f'(z)+df(z) = 0, /(0) = 0, f(0)= 7 . (22) 



Denoting g — f one can rewrite Eq. ( f22"|) as a first order differential equation for 
9(f)- 

g^ = f[2g(f)-d], g(f)\ f=0 = r (23) 



If 7 < d/2 solutions to problem (|23|) are given by the following periodic trajectories 
in the phase space (/, /') = (f,g): 



, d d/2-g 



f 2 =g-l + -\n 



2 \d/2-jj' 

In particular, for \ j\ d the trajectories are given approximately by ellipses 

d-f + g 2 = 7 2 . 

This family of regular solutions corresponds to the straight line r] = 2 — d for all 
7 < d/2 on the (7, ?7)-plane. Note that for d = 2 they are solutions with rj = 0. 

If 7 > d/2 the trajectories in the phase space are unbounded. They describe 
solutions f(z) which are singular at some finite value of z. 

Let us consider now the case A + = 0, i.e. r\ = 2 + d. Eqs. (|20| ) and (^1|) become 



f"(z)-2f(z)f'(z)-dzf'(z)=0, /(0) = 0, HO) =7. (24) 

Solutions of this problem can be obtained from solutions of Eq. (^) with 7 substituted 
by 7 + d/2. Indeed, let function h(z) be a solution of 

ti'(z) -2h{z)h'{z) +dh{z) = 0, /i(0) = 0, /i'(0) = 7 + - 

2 

(cf. (|2~0"D). As one can easily check, then f(z) = h(z) — {d/2)z is a solution of problem 
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In a general case solutions of (|20|) , (|2T|) can be found only numerically. A generic 
solution f(z) or its derivative end up with a singularity at a finite point z = z , i.e., 
is not regular. A natural method for searching regular solutions is to fine tune rj and 
7 in such a way that the position of singularity zq = zq^,^) — > oo. A realization 
of this recipe in a practical computation is to pick a solution for which zq shows a 
tendency to diverge |2^, j2~6fl . 

First we studied numerically solutions of the FP problem, Eqs. ( PDj ) and for 
d — 2, in a wide range of values of 7, r/. We found that for rj = and any 7 < 1 
there is a periodic solution. This is in accordance with our analytic result for A - = 
above (note that for d = 2 A~ = — f]/2). 




Figure 1: Non-trivial FP curves for d = 2 and d = 3. Only the first 7 curves are 
shown. 

For < rj < 1 and rj < 2 — 27 the parameter 7 can be fine tuned to a certain 
value so that a non-trivial regular FP solution exists. The set of such values (7, if) 
form an infinite discrete series of curves 77^(7) (n = 1, 2, 3, . . .) shown in Fig. |I[ At 
7 = r) n (0) = 2/(n + 1). The curves accumulate at the line r\ = as n — ► 00. 

There are general grounds to expect the existence of the curves r/ n (7) whose points 
(7, //n(7)) correspond to regular FP solutions p2| . Indeed, let us suppose that for some 



value (7', rj') there exists a regular solution. This means that the singularity of this 
solution is situated at ^0(7'; v') = 00 • Consider now another value rj" sufficiently close 
to rj'. Assuming that the function zq^,-)]) is continuous in the vicinity of (7', 7/) it 
is clear that there exists the value 7" such that again ^0(7", 77") = 00. Consequently, 
there is a curve 77(7) passing through the point (7', r/) in the parameter space. The 
curve is defined by the equation 2:0(7, v) = °°- 

We found that when moving along a given curve f] n {l) the shape of the corre- 
sponding solution qualitatively remains the same. When passing from one curve to 
another the shape of the solution changes significantly following the same regular 



pattern as the one obtained by T.R. Morris within the Legendre ERG equation [E7j. 
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Since all the curves 77^(7) are situated in the rj > half-plane, their presence is a 
signal of the existence of an infinite number of non-trivial FPs in two-dimensional 
scalar theories. This expectation turns out to be true. As it was shown in Refs. P7[ 



and |32| by studying the next-to-leading ERG equations, they do correspond to the 
minimal unitary series of p(p + 1) (jp — 3, 4, . . .) conformal models. 

For d = 3 we found periodic solutions corresponding to 77 = — 1. Again, this is 
in accordance with our analytic results for A~ = because for d = 3 the condition 
A~ = leads to 77 = —1. We also found an infinite number of critical curves rj = 77,1(7), 
n = 1, 2, . . ., corresponding to non-trivial regular FP solutions of Eqs. ( PU[ ) and (|2~T|). 
They are plotted in Fig. |l[ The curves accumulate at the line 77 = — 1 asn- * 00. At 
7 = rj n (0) = (2 — n)/(n + 1). The curve 771(7) crosses the 7-axis at 7* = —0.229 . . .. 
The value of 7 coincides with the one for which a non-trivial physical FP was found 
in the LPA of the Polchinski equation fl26| . This is the well-known Wilson-Fischer 
FP analyzed in numerous previous studies (see P,[p2]]-[p^], plf , |3^]-[0. It belongs 
to the Ising model universality class. The rest of the curves are situated between the 
horizontal lines 77 = and rj = —1 (except for the value of 1/2(7) at 7 = 0: 1]2(0) = 0, 
the GFP). This suggests that for d = 3 there is only one non-trivial physical FP, 
namely the one associated with the curve 771(7). 



3 N = 00 Polchinski ERG equation 

In this section we will study solutions of the leading order Polchinski ERG equation 
for N = 00. In this limit many features of the solutions can be described analytically, 
and we are going to take advantage of this. Our aim is to study the 77(7) curves of 
regular FP solutions for N = 00. We expect them to be certain smooth deformations 
of the analogous curves for finite N, therefore the analysis of the N = 00 case will 
help us to understand the features of the 77(7) curves for finite N. 



3.1 Equation and general solution 

In the limit iV — > 00 the Polchinski flow equation, Eq. 
and becomes 



5D, simplifies considerably 



ii(y, t) = -2yuu' - u 2 + ( 1 + 2yA 



u y + su 



(25) 



for —00 < t < 00 and y > (see Eq. (|T3|)). 

The FP equation for u(y) follows from Eqs. (17) and (|T8|), and is given by 



2yuu' y - u 2 + (l + 2t/A~) u' y + su 

u(y)\y=o 



0. 

27. 



(26) 
(27) 
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As in the case of finite N, Eq. (|26|) is solved by the GFP u{y) — and by the TFP 
u(y) = s. 

Let us consider solutions of flow equation (|2lS| ) first. As was pointed out in Ref. pSfl , 
it can be solved analytically for the inverse function y(u,t). Indeed, from Eq. ( p5|) 
one obtains the following flow equation for y(u,t): 

y{u, t) = {u- s)uy' u - 2 (A" - it) y - 1. (28) 
Its general solution is equal to |^9| 

!/(«. *) = 7—^37^ f 1 ^^ + y(ti), (29) 

where 



[s - u) 2 - a u a ' \ u 



a = 2A~/s, (30) 

F is an arbitrary function and y(u) satisfies the inverse FP equation, i.e. Eq. ([28|) 
with y — 0, and the corresponding initial condition: 



(tt-s)Vu-2(A--«)3/(«)-l = 0, (31) 

2/(27) = 0. (32) 

The FP equation, of course, can also be obtained directly by inverting Eq. fl26|) . 
Condition follows from (g?p. 

Eq. (|3T|) can be solved analytically p9| . Its general solution is given by 



VW 



1 



u) 2 - a u a 



u a {s-u) l - a a - 1 r (s - ^\- Q , 



(33) 



Note that it is consistent with general solution ( p9]) with the function F taken to be 
constant, F = —C. 

The form of the solution is not unique and integrating by parts one can bring it to 
another one. Possible divergences of the integral in Eq. (|33l) can be absorbed in the 
integration constant C . Equivalently, one can choose the lower limit to be uq 7^ 0, in 
this case C becomes uo-dependent. 

It is important to bear in mind that a solution u(y) may not be monotonous and, 
therefore, not invertible on the whole semi-axis y > 0. As an example let us consider 
the case in which u(y) is monotonically decreasing in a subinterval < y < y# and 
monotonically increasing in the infinite interval y# < y < 00. One has u(0) > u(y#) 
and u(oo) > u(y#). In each of these regions u(y) is invertible, let us denote the 
corresponding inverse functions as y^(u) and y(>)(u) respectively. The former is 
defined for u# < u < u(0), where w# = u(y#), the latter for u# < u < u(oo). They 
have to satisfy the following matching condition: 
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y {<) {u # ) =y {>) {u # ). (34) 

Hence, there are two branches of y(u), y^(u) and y^(u), such that their inverse 
functions form the solution u(y) for < y < oo. Note that initial conditions ( |2"7| ) and 
(|32| ) imply that ?/(<)(m(0)) = y(<)(27) = 0. This equality fixes the integration constant 
C in formula Q3"3] ) for the ?/(<) (w) branch. The analogous constant in the expression for 
the ?/(>)(«) branch is fixed by the matching condition, Eq. ([34]) . Concrete examples 
will be considered in Sect. |3~3| 

To conclude this discussion we stress once more that Eq. (|33| ) gives the exact 
general solution of the Polchinski FP equation, Eq. (|26|) , in the LPA in terms of the 
inverse function y(u). The constant C is fixed either by initial condition fl3"2|) or by a 
matching condition similar to ([34]) . 

The parameter a, defined by fl3"U|), turns out to be very convenient in the analysis 
of FP solutions. The anomalous dimension rj and the parameter s are expressed in 
terms of a as follows: 

d d . . 

77 = 2-- , s =- . 35 

1 — a 1 — a 

There are different types of solutions for different intervals of values of a. The interval 
which is most interesting for physical applications is 

- 1 < a < 0. (36) 

It corresponds to the range 2 — d < 7] < (4 — d)/2. From now on we will only be 
considering interval ( |36D of values of a. The form of the solution, given by Eq. (|33|) , 
is adequate to this interval. 

The asymptotics of the solution u(y) for large y is given by the N — > oo limit of 
Eq. (|19D. In the case in hand it can also be found from exact formula (0). For this 
we first calculate the asymptotic behavior of the function y(u) as u — * s and then 
invert the asymptotic formula. This procedure not only reproduces formula (|19|) , but 
also allows to determine the constant D. One obtains 

D 



Sin 7T( 

Now we turn to the analysis of the expression for the derivative of the solution 
u(y) following from Eq. (f26[): 



u(s — u 
2uy - (1 + 2yA 

The denominator is zero if 2yu = (1 + 2yA~). Let us introduce the function 



U y = — M , n-A-V ( 3T ) 



1 + 2yA~ 

U s ing(y) = • (3£ 
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For a satisfying ( |36]) we have s > and A~ < 0. It follows that a general solution u(y) 
with initial condition u(0) = 27 (see fl2~7|) ) and asymptotics lim^oo u(y) = s (from 
Eq. (|T§|)) necessarily crosses the curve of potential singularities u sing (y) at some point 
y = y S in g - In other words, there always exists at least one point y S ing > such that 
u{y S ing) = u sing (y sing ). Let us denote this value of u as u Q , 

1 \ 1 ~l" ^y sinq^- 

= u(y sing ) = 



2y. 



sing 



Since in general « 7^ 0, «o 7^ ^ the derivative of the solution does not exist at this 
point and the solution is singular .0 

Note that this picture has an equivalent description in terms of the inverse function 
y{u). Its derivative is obtained from Eq. ([31]), 



, _ l + 2(A~-% 
Vu u{u - s) ' 

So the line of potential singularities is 

Vsingiu) = 2(M _ X A _ r (39) 

the inverse function of u sing {y). 

A general solution y(it) always intersects this curve at some point, namely at 
u = uq (see Fig. |2|). Note that in general there may be a few points of intersection. 
If u 7^ 0, u 7^ s the derivative y' u vanishes, the derivative u' y does not exist, and the 
solution is singular. 

The analysis of relation fl37|) shows that a general solution has the following be- 
havior near the point of singularity: 



u(y) =u + A^Jy- y sing + b(y - y sing ) + c(y - y sing f /2 H . (40) 

The coefficients are determined from Eq. (|37| ) and are equal to 



A = =h 



U (s— Uq) S + A — 3u 

etc. 



\ ysing 3?/ 



sing 



Using these expressions one can easily see that the solution is indeed non-analytic 
and the derivative does not exist at y = y S i ng : 

Av) = , = = + b+---. (41) 

^yV ysing 

By inverting a general solution, Eq. fl4"0|) , one gets the following expansion for y(u) in 
the vicinity of u — u : 

2 This observation was made in Rcf. |29|. 
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Figure 2: Plot of the curve of potential singularities y S in g (u) defined by Eq. (|39|) for 
d = 3 and A - = —3/8. Also given, respectively as bold and dashed curves, are the 
plot of a singular solution y(u) satisfying the initial condition 2/(27) = 0,7 = —5/2 
and the plot of a regular solution y(u) satisfying 2/(27*) = 0,7* = —1/2. The point 
where all these curves intersect is (0, y if = —\/ (as)) with a = — 1/3 and s = 9/4. 



vw 



ysing ~l" ^2 



2k 

\ 2 / \ 3 

u o) -^( u - u o) + 




U-Uq) + 



The corresponding derivative is 
2 , v 6S 



A 2 



it 



M0 J 



A 4 



mo) 2 + 4 




M ) 3 + 



which at u = u is indeed zero. 

The situation changes if the parameter 7 in the initial condition, Eq. (0), is 
adjusted to a value 7* such that u(y) crosses the curve of potential singularities at 



y = y*, where is determined by the condition u(y* 
one can see that 



u 



sing 



(y*) = 0. From Eq. 



y* = - 



1 

as 



(42) 



Such a solution u(y; 7*) does not have root behavior and u'(y) is regular at the point 
of potential singularity y = y*. Consequently y' u is no longer zero at u = (see Fig. ^| 
where a plot of y(u; 7*) satisfying j/(0; 7*) = is also given). 

The expansion of a regular solution in the vicinity of y = contains only positive 
integer powers of (y — y*) with the leading term being 
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u(y) = b n (y - y*) n + higher powers of (y - y*). (43) 

It is the leading power n which defines the type of the solution. It turns out that 
there are two qualitatively different classes of regular solutions: (1) with n = 1, and 
(2) with n > 2. 

For n = 1 the expansion at y = y„ is given by 



u(y) = b 1 (y-y*)+b 2 (y-y*y + b 3 (y-y*y + ---, (44) 
s 2 a(a + l) 3s 3 a 2 (a + l) 2 

61 = T- 1 ' h = ~ 4(a +2) - (45) 

For n > 2, as one can find out using ERG FP equation (|26|) or Eq. fl37|) , the behavior 
of the solution is more complicated: 



u(y) = bniy-y^ + c^y-y^ + hniy-y*) 2 " (46) 
+ d 3n (y - y*) 3n ~ 2 + c 3n (y - y^ 1 + ■■■. 

For such a function to satisfy the equation the parameter a must be fixed to the value 
a = —1/n. This amounts to 77 = 2 — dn/(n + 1), s = dn/(n + 1), etc. However, the 
coefficient b n is not fixed by the equation and remains a free parameter. In other 
words, there is a regular solution for any value of b n . The rest of the coefficients are 
expressed in terms of b n . For example, 

C 2 n = Tl 7T & «' &2 « = & n; etc - 

s z [n — 1) s 

Note that for n = 2 the terms proportional to (y — y*) 2n and to (y — y*) 3n_2 give the 
same power of {y — y*) 4 so that expansion (|4*6| ) is of the form 

u(y) = b 2 (y - y*) 2 + c 4 (y - y*) 3 + (64 + d e )(y - y*) A H . 

Let us analyze the expansion of the corresponding inverse functions ?/(m) around 
y*. Consider the two classes of regular solutions. 

(i) For the solution u(y) with n = 1 the inverse function y(u) has the following 
behavior close to u = 0: 

y{u) = + l u - | M 2 + 0( M 3 ), (47) 
where 61, 6 2 are given by Eq. (f4~5|). 
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(ii) For regular solutions with n > 2, = n/s and the expansion for y(u) is the 
sum of two series in powers on u of the form 



Tl i r i 

y(u) = — h cm™ 1 + k x u + k 2 u 2 + . . . + l\u + l 2 u 2 + • ■ ■ , (4£ 
s L J 

where a, k iy lj are related to the coefficients in Eq. (|46| ) as follows: 

1 6 2 « 2n + 1 c 2n 2n 2 

k\ = — = , h = — = —. — , etc. 



l)H n '' nb\ sn ' nb\ s 2 {n — 1) ' 

This analysis explains the nature of the condition of regularity of solutions of the 
ERG equation (for iV = oo): the parameter 7 must be adjusted to the value for which 
the solution u(y), fixed by the initial condition u(0) = 27, satisfies 

u (-—J—) = 0. (49) 



2A" 

This property actually follows from Eq. (pT|). Indeed, it is easy to see that if the 
derivative y'{u) is finite at u — 0, then y(0) must be 



The latter is equivalent to ([49]). 

We would like to note that there is another class of solutions, namely those which 
satisfy 

u{yl) = Usingiy'J = s, 

where y^ = 1/(2A + ). As it can be seen from (|37D and (^), its derivative does not 
have a singularity and the solution is analytic at the point of potential singularity. 
This class will not be discussed further in the article. 

From now on we focus on regular solutions of Eqs. (|26|) and (p7|) only. When 
it is necessary to indicate their dependence on the parameters we will be using the 
extended notation u(y; 7, r], d) and y(u; 7, i], d) correspondingly. 

Expanding formula (B3f) in powers of u around u = one gets 



y[U) = 

as 



C 1 

s 2-a u a 



1 - 



a - 2 



u + 0{u 2 



a(l + a)s' 



: u + 0{u 2 ). (50) 



Compare this expansion with Eqs. (|47D and fl48|) . It is easy to see that there are two 
possibilities for function (|50"D to give a regular FP solution: 



(i) C = 0, in which case solution (pE?D or (f44"|) is reproduced 
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(ii) (7 / 0, a = — 1/n, n = 2,3, . . ., in which case solution (f48|) or (f46|) is obtained. 



One may also wonder whether there is a regular solution in the case C / 0, 
a = — 1. It can be readily checked that the answer is negative. Indeed, for a = — 1 
the integral in (|33|) can be calculated explicitly. To avoid the divergence at the lower 
limit let us consider a modification of Eq. (^) with the lower limit of integration 
u = substituted by u = u > and the constant C changed to another constant 
C uo . Performing the calculation one gets 

, . u 1 2u(u — uo) 2siu u . ri . 

(si — u) 6 si — u [si — u) 6 {si — uy u 

where S\ = d/2. For a = —1 it follows from Eq. (^5|) that rj = (4— d)/2. The behavior 
of solution fl5T|) at u = is given by 

1 2 

y(u) = 1 — kU \nu + 0(u,u 2 Inu, . . .). 

St s( 

Consequently though u{l/s\) = u'{l/s\) = the second derivative u"(y) is singular 
at y = l/s\ and so the solution u(y) for a = — 1 is not regular. 

In the following subsections we will see that the two classes of functions listed 
above are quite different and characterized by two distinct sets of eigenvalues and 
eigenoperators for perturbations around the FP [25|. Let us consider both cases in 
turn. 



3.2 C = 

In this case the solution u(y) is invertible for all < y < oo. The integration constant 
in Eq. ( |33D is fixed by initial condition ( |32"D and, therefore, for a given d becomes the 
following function of 7 and rj: 

a a Jo \ z J 

The condition C(j, 77) = defines a curve in the parameter space which we denote 
as 771(7). The lower index "1" corresponds to the power n — 1 of the leading term in 
expansion ([43"D. Using Eq. (|52|) the condition (7(7,77) = can be written as 

- ^(1 - a) 2 (l - ^(1 -a))" 1 f 1 dzz a (l - ^(1 - a)z) " = 1, (53) 



d V d J Jo \ d 



where we took into account the second relation in (p5]) . Solutions in the case un- 
der consideration will be denoted as Ui(y;j, d) = u(y; 7, 771(7), d) and yi(u;j, d) = 
y(w; 7, 771(7), d). They exist only for negative values of the parameter 7. This can 
be seen from the following arguments. Eqs. (|35|), (|36|) and (|42| ) imply that s > 
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and > 0. In accordance with asymptotic formula ( |19D the solution u\ (y;^,d) ap- 
proaches s for large y. As it follows from Eqs. (|44|) and (|45|) , Wi(t/*; 7, d) = and has 
a positive linear slope at y = y*. Therefore, regular solutions exist for 7 < only, 
since otherwise there would be other zeros of ui(y) in the interval < y < and 
u i(y*, 7; d) would not be invertible. 

Thus, Eq. ( |53"D defines a as a function of 7 for 7 < 0. In fact, one can easily see 
that it is a function of the ratio j/d. We denote it by 0:1(7/6?). In accordance with 
the first relation in Eq. ([35|) this defines the function 

»/i(7) = 2-- ^y-v- (54) 

Limiting values of 0:1(7/6?) can be readily derived from Eq. (^). For 7 < and 
7^0- 

«i [ 2 ) - -1, or 771(7) - — . 

These limiting points correspond to the GFP. Expanding Eq. (|53|) in powers of 7 and 
using Eq. fl5"4] ) we obtain that 

m (7) = — + 2 7 + 24^- + 480^ + d ■ O Uj . (55) 



For 7 — ■* —00 



Q) ^0, or 7^(7) ^2-d. 



The asymptotic expansion of the function 01(7/6?) can be found from the analysis of 
Eq. (|53D. Introducing the variable w = —6?/ (27), which is positive and tends to zero 
as 7 — ► —00, we get 

Oi ( ] = — w + w 2 — w 2 In w — w 3 In 2 w + • • • . 



The dots stand for the terms w 3 , w 4 , w A \n 3 w, etc., which in this limit are at least 
logarithmically smaller than the terms written down explicitly. 

As before, we denote by 7* the point at which 771(7) = 0. As was explained in the 
Introduction, this point is of special interest because it corresponds to the physical 
FP solution of the Polchinski ERG equation in the LPA. From expansion (|55D one 
sees that 7* = for d = 4. Hence, in this case only the trivial GFP exists. For d = 2 
we find 7* = —00 and, therefore, there are no non-trivial physical FP solutions for 
finite values of 7. 

Let us analyze the case d = 3. The value 77 = corresponds to o = —1/2 (see 
Eq. (0)). The integral in formula can be calculated explicitly. One gets the 
equation 
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-v / 2^(2 - 2k,)- 3/2 (sinh(20*) 
where the parameter 9* is defined by the relation 

/ K$ 

~2 ' 

and = e~ llT 2^* > 0. 

d=3 



29*. 



(56) 



sinh 9* 



Figure 3: Plot of the function 771(7) for d = 3. 



To obtain the curve 77 = 77! (7) we solved Eq. fl53|) numerically for 7 < 0. For d = 3 
the result is given in Fig. ^. In other dimensions the function 771 (7) has a similar profile 
and can be easily computed from the d = 3 curve using the scaling properties which 
will be obtained in Sect. |j. The physical FP solution u*(y) = u\(y; 7*, 3) corresponds 
to 7* = —0.634913 . . . which is the value satisfying Eq. ( |56|) . The function u*(y) and 
its inverse y*(u) are given in Fig. f|. The numerical results for the curve 771(7) confirm 
that non-trivial regular solutions exist for 7 < and — 1 < t] < 1/2 only. For 7 = 0, 
77 = 1/2 we have the GFP. Recall that the same is true for N — 1. Note as well that 
the curves 771(7) for N — 1 and N = 00 have very similar shape (see Ref. |32[] and 
Sect. H). Since the curve 771(7) has a part with positive values 77 > it corresponds 
to a physical FP. 

The eigenoperators of perturbations around the FP solution u = 161(2/57, 3) are 
equal to 



9m(y) = ex P 



where 



dz t r 



(57) 



A, 



A m + 2yu' y — s + 2u 

1 + 2A~y - 2yu 
—s(m + a). 



(58) 
(59) 
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d=3, eta=0 d=3, eta=0 
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Figure 4: Plots of the functions Ui(y) and y\(u) for d = 3 and 77(7*) = 0, 7* 
-0.634913 .... 



The spectrum (|59| ) and the linear behavior of u*(y) at y — y*, Eqs 
indicate that the physical FP at 7 = 7* is the Wilson-Fischer FP p9[ 



and 



3.3 C^0, a = -l/n 

Now we consider regular FP solutions with C 7^ and 



n = 2,3, 



n 



in Eq. (50). From relations fl35|) we get that in this case 



dn 

n + r 



1 — a. 



dn 
n + 1' 



(60) 



In what follows we will use the notations u n (y; r y,d) and y n (u;j,d) for the FP 
solutions u(y;j,T] n , d) and y(u;j,T] n ,d) respectively. It follows from Eqs. ( fIBl) and 
( |50| ) that the leading behavior of the solutions in the vicinity of w = is given by 



y n (u;7,d) = y m + 



C 



2+1/n 

Sn 



and in the vicinity of the point of potential singularity y = y m = n/s n = n{n + l)/d 
is equal to 

g -(2n+l) 

u n {y, 7, d) = n Cn {y - y m ) n + ■■■. 

The curves of regular solutions are straight lines 77 = r] n = const in the (7, r?)-plane. 
We consider the cases of even n and odd n separately. 
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Figure 5: Plot of the function 6*2(7) for = 3. 



3.3.1 n even 



In this case a solution y n {u\ r y ) d) has two branches. As it was discussed earlier, 
the branch which satisfies < ?/„(«; 7, <i) < y m is given by formula ([33|) with the 
integration constant equal to 



CM 



(2 1 ) a -(s n -2 1 ) 1 - a - a n -l r*yfs n -z 



a, 



OC n JO 



dz 



(cf. (|5"2D). This expression is a consequence of the initial condition, Eq. 

For n even the function C n (j) and the solution y n (u; 7, d) exist for < 7 < s n /2. 
In particular, for n = 2 one gets 



C 2 {l) = +3s 2 arcsin J— + 2a j 2 7 (g 2 + 7), 

with S2 = 2d/3. The plot of the function 6*2(7) for d = 3 is presented in Fig. |5]. 

For each value of 7 > and 7 < s n /2 the two branches of the solution y n (u; 7, d) 
are given by: 

for < u < 27, < y n (u) < y* n 



vw 



CM 



1 a - 1 1 

+ 



s — z 



u a (s — u) 2 a a(s — u) a u a (s — u) 2 a Jo \ z 

u 1/n r( \ n (n + l)u 1/n [ u (s n -z\ 

~{s- M )2+i/™ ~ ( Sn _ u) + {s n - M ) 2 +V« J Q 

for < u < s, y^n < y n {u) < 00 



dz (61) 



dz: 



cm 



u a (s-u) 2 - a a(s-u) 



+ 



a — 1 



a u a (s - u) 2 ~ a jo 



s — z 



dz (62) 
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u 1/n n (n + l)u l/n /•" fs, - 



u 



Cud) -j— t+ , _' /n I ^— dz. 



] 2 +V" HKU (S n -U) (S n - U) 2+I l n Jo V Z 



Formulas (KTj), (p2"D match together at u — 0. These two branches define a unique 
function w n (y; 7, d) for 7/ > 0. In Fig. |] we show the plots of the functions 7x2(7/; 1/2, 3) 
and y2(u; 1/2,3) respectively. 

It is easy to see that C n {^) — > +00 when 7 — > + . This limit, of course, corre- 
sponds to the GFP. For 7 = s n /2 we find 

6 n (s n /2) = nd— 



sin to„ sin 



The corresponding regular solution is u n (y;s n /2,d) = s n = const, the TFP. This 
function is not invertible, hence the corresponding solution y n (u;s n /2,d) does not 
exist at 7 = s n /2. 



3.3.2 n odd 

For n odd, n > 3, regular solutions exist only for negative 7 and consist of one 
branch. The function y n (u\^,d) is defined for 27 < u < s and varies in the range 
< y n (u; 7, d) < 00. The integration constant C in formula ( |33D is fixed by the initial 
condition, Eq. (|32|) , and is given by Eq. (|52|). As before, we denote it by C n {j). It is 
easy to see that for odd n as k = e -i7m 27 — > + 

C„(7) ~ -»■ -00, 

The plot of the function 6*3(7) for c? = 3 is shown in Fig. |7[ The solutions u 3 (y; —1, 3) 
and yz{u] —1, 3) are given in Fig. |[ 

Note the characteristic cubic root behavior of 7/3(7/; —1, 3) in the vicinity of u = 
(see Eq. @). 
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Figure 9: Plot of the function ys (u) for d = 3 and 7 = —2.5. The function is 
non- monotonous in a vicinity of u = 0. 

At certain value 7 = 7* n < the function C n (j) =0. Of course, this is the 
point where the straight line rj — r\ n (n = 3, 5, . . .) in the (7, r/)-plane hits the curve 
V — Viil) obtained in Sect. |3~T| . 

The solution u n (y; , y* n ,d) coincides with ui(y;j* n , d) for this value of 7. When 
trying to search for solutions for 7 < 7* n one finds out that y n (u; 7, d) becomes 
non-monotonous, and hence non- invert ible, in a vicinity of u = 0, and therefore the 
solution u n (y; 7, d) does not exist. An example of such behavior is illustrated in Fig. || 
which shows the plot of the function ys(u; —2.5, 3). 

The eigenoperators are of the same form as in Eqs. fl57|) and (|58|) . The eigenvalues 
are given by 

Tfi 

A"= Sn (l--). (63) 
n 

Note that these are also the eigenvalues of the GFP. 

Let us summarize the results obtained in this section. We have found the infinite 
discrete series of lines r\ = rj n {j) {n = 1,2,3, . . .) located in a certain region of the 
parameter space (7,^). Points on these lines correspond to regular FP solutions of 
the N = 00 Polchinski ERG equation in leading approximation. The curve 771(7) (see 
Fig. ^) bounds the region from the left. From the right the region is bounded by the 
straight line rj = 2 — 27 of TFP solutions u = s. The curves rj = rj n , n = 2, 3, . . ., 
are straight horizontal lines (see Eq. (0)). They start at 7 = (GFP solution) and 
end up either at 7 = s n /2 > (the TFP solution) for n even, or at 7 = 7* n < for 
n odd. The value 7* n is defined by the condition r7i(7* n ) = Vn- As n — >• 00 the lines 
f] = Tj n accumulate at r/oo = 2 — d. 
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4 Solutions for finite N and their N —> oo limit 

In the sections above we have seen a detailed treatment of the Polchinski ERG equa- 
tion for N = 1 and N = oo. Not presented there are the finite N cases, that are of 
interest in their own right, but also because the N = oo case should be viewed as a 
limiting case for large N. For that reason, we present in this section a study of the 
finite N case, and show how the N = oo case comes about as a limit. 

First of all from the plots for N = 1 in Fig. [I] it is easy to notice that the patterns 
of the curves of regular solutions are the same for d = 2 and d — 3. They differ just 
by a vertical shift. Our analysis in the previous section shows that this is also true 
in the N = oo case. It turns out that for a given N the pattern of the curves 7^(7) 
is universal and does not depend on the dimension of the space d. Namely, when we 
pass from one number of dimensions d to another d the functions 77^(7) experience 
a constant shift and some scaling transformation while the pattern of the curves is 
preserved. This is a consequence of a property of the ERG equation which we are 
going to analyze in the remainder of this section. 

Let u(y;j, 77, d) be a regular solution of fll7D , fll8D corresponding to a point of a 
curve 77 = f] n {l)- Perform a scaling transformation 

u(y) -*u(y) = \u(y;r),j,d), V=\v- 
It is easy to check that the function u(y) satisfies the equation 



2y d n ~~du o / 2 „„. A du . _, „ ,„ jS 

iw U(V) ~ 2VU d§ ~ \ ~N ~^ 2yA~AJ — + sXu — 0, (64) 

u(0) = 27 = 2A7. (65) 

Therefore u(y) is a regular solution in d dimensions corresponding to the point (7,77) 
of the curve fj = fj^j), where 7, 77, d and a relation between rj n {^) and 77 n (7) can be 
found from Eqs. (Q) and (p5]) as follows. According to Eqs. (^) and QT5|), r\ = 2 — s, 
d = s — 2A~. Correspondingly, fj = 2 — s, d = s — 2A _ , where A - = AA~, s = As, 
see Eq. (0). Combining these relations we obtain that 



d = \d, 7 = A7, 

fj = 2- As = 2- A(2-?7) = 2(1 — A) + A77. 
Trading A for d/d we arrive at the relation 

This formula is valid for any N > 1. 
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Eq. (|66|) can be rewritten as the relation 



fjjjij) - 2 _ rjnjdrf) - 2 
d ~ d 

which tells that the combination (3 n {^) = (r] n (dj) — 2)/d is in fact independent of d. 
Therefore, a general solution of functional equation fl66[) is given by 



^( 7 ) = 2 + d/3 n Q). (67) 

The functions 77^(7) studied in the previous sections do have this property. We 
checked that in the case N = 1 and d = 2 or d = 3, the curves 77,1(7) verify relation 
( |5"6j ) or, equivalently, are described by flB7|). Special limits A + = 0, A - = are 
consistent with ( |66|) , (|67|) as well. The results for the N 00 case in Sect. 3 are also 
in a full correspondence with this property. In particular, from Eqs. (54) and (SO) it 
follows that 



0n 



7 



dJ ~ 



Recall that for n > 2 a n is the constant function a n = —1/n. 

For d fixed the shape of the curves changes with N. For N = 1, d — 2, 3, a few 
curves are shown in Fig. [I]. For iV = 00, d = 3 the plot of 771(7) is given in Fig. [3[ 
For n > 2 T] n = 2 — dn/ (n + 1) = const, i.e., they are represented by horizontal lines. 
It is of interest to study the curves for finite N > 1 and understand the transition 
from iV = 1 to N = 00. In particular, a natural question arises: Do the intersections 
of the curves 7^(7) (n odd, n > 3) with the curve 771(7), encountered for iV = 00, 
have an equivalent for finite iV? Such intersections do not exist for N — 1. Is there 
a value of N for which intersections are encountered, or the N = 00 curves approach 
the large iV limit in some other way? 

To gain more insight into the properties of the functions 77,^7) for an arbitrary 
N it is useful to study their behavior in the vicinity of 7 = 0. This can be done 
by a perturbative analysis for small values of 7. To this end, perform the change of 
variables v = —NA~y, u(y) = 2jh(v). Eqs. fllTD and ( |18D turn into 



L(v)h-Eh + ^vhh' z + -^h 2 = 0, (68) 
h(0) = 1, (69) 

where 

E= 2~K Z ' L{v)h = vh' m + {M-z)ti v 

and 
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The linearized equation L(v)h = Eh can be solved using the known solutions of 
the eigenvalue problem L{v)h n = X n h n : 



X n = -n, h n (v) = $(-n,M;v), n = 0,1,2, (70) 



where 



-r / 7 \ 1 a (0)2 v 2 

0(o , M = 1 + r + LL_ + ..., 

(o) n = a(a + 1) . . . (a + n — 1), (a) = 1, 

denotes the confluent hypergeometric function (see Ref. [fig] ). Note that /i n (t>) is a 
polynomial of order n with /i n (0) = 1. 

We now develop perturbation theory around these solutions. Fixing n, we write 



h(v) = (H 7 D 1 + 7 2 I) 2 + ...)(/ i » + 7E c A + 7 2 E (| A + '' 



(here V n are constants) and solve system (|68|), (|69^) order by order in 7. Note that, 
at any order, the sums over k 7^ n are over a finite number of terms. Here we list a 
few results. 

At order 7 , we find E n = E^ = —n, from which it follows that 

t\i d , , , , 

^ = 2-^ + 0(7), (71) 

a result that is independent of N. A more detailed calculation yields 



770(7) = 2-2 7 , 

d , 3 . ,24, 13 12, . o 

*W = 2 -2 + 2 ^ + M> +7 V 1+ M + AP )+0h 
2d 3M + 8 /rw 2 , 

= T ~ 7 m(mTT) + * 7 

3d „ 3M 2 + 69M + 196 ,^ 2 , 
%(7) = 2 -T + 67 M(M + l)(M + 2) +0(7) ' 

, , 4d 15M 2 + 215M + 636 ^, 2n 

^ = If" 487 M(M + 1)(M + 2)(M + 3) + 0(7 } ' 
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776(7) 



5d 



Qd 

2 

7 



5M 3 + 360M 2 + 4225M + 12846 2 
M(M + 1)(M + 2)(M + 3)(M + 4) + ^ 7 

35M 3 + 1365M 2 + 14308M + 44364 2 
96 ° 7 M(M + 1)(M + 2)(M + 3)(M + 4)(M + 5) + ^ 7 



Note that the expression for 770(7) is, in fact, exact as presented, i.e. it does not 
have 7 2 corrections. One recognizes it as the line of the TFP solutions. Moreover, 
note that the order 7 contribution to r/ n (7) vanishes in the limit M — > 00 (and thus 
for iV — > 00) for n > 2. This is, of course, consistent with what we expect, namely, 
r] n constant for N = 00 and n > 2. We would like to stress that the expansions for 
Vn{l) obtained here are in a full agreement with general formula fl5"?D expressing the 
universality of the pattern of the curves for various d. 

To obtain the curves 7^(7) for arbitrary values of 7 a numerical approach was 
taken. We calculated them for a wide range of values of N (up to 1000). As an 
illustration let us consider the case d — 3 in more detail. For even n, the curves lie 
in the region delineated by the lines 7 = 0, rj = — 1 (7 > 0), and rj = 2 — 27 (this is 
the curve 770(7))- For odd n the delineating lines are 7 = 0, rj = — 1 (7 < 0). Higher 
values of n correspond to lower-lying curves. For 1 < N < +00 the pattern of the 
curves is very similar to the the N — 1 case, plotted in Fig. |l|. 

We found a one-to-one correspondence between the curves for various iV and fixed 
n. In Fig. [10] we have plotted the case n = 4, where the 770(7) line is also included. 
We see that for increasing values of N the curves lie higher. In the N —>■ 00 limit for 
< 7 < s n /2 they approach the (horizontal) line rj = rj n and for s n /2 < 7 < 3/2 the 
line 770(7). Similar curves but with lower precision were found for odd n > 3. 



-0 . 8 




Figure 10: Curves in the (7, 77)-plane representing regular solutions for the case d = 3, 
n — 4, with iV = 1, 2, 5, 10, 20, 50. Curves lying higher correspond to larger values 
of N. The dashed line is rj = 2 — 27. 
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Within the accuracy of the computation it can be concluded that as N grows the 
functions r) n (j), existing for — oo < 7 < 0, transform to r] n =const, existing only for 
7*n < 7 < 0. The latter is the limit we expect from our studies of the case N = 00 in 
Sect. 3. 



5 Discussion and conclusions 

We studied families of regular FP solutions of the Polchinski ERG equation for the 
0(iV)-model in the LPA. The families are labeled by the integer n > 1 and are 
represented by curves ?/ n (7) in the (7, ?/)-plane. We proved that for given N the 
pattern of the curves is universal for all d and described by Eq. (JJTp. 

There is a one-to-one correspondence between the curves for different N. In 
particular, for any N 

Vn(0) = 2- '"' 



n + 1 



see Eq.([n]), and rj n — > 2 — d as 7 — ► —00 for n odd. Further properties of the curves 
were discussed in detail in Sect. f|. 

As N increases the 1/1(7) curve transforms into the corresponding curve for N = 
00. The rest of the curves transform into the straight lines given by Eq. (pU|). 

For the N = 00 case the solutions were studied analytically. We have shown 
that most of them have a singularity at some finite point (see Eq. (pH])) and are not 
acceptable from the physical point of view. The condition of regularity imposes a 
relation between the parameters rj and 7. We analyzed this condition of regularity 
explicitly and demonstrated how various classes of FP solutions appear. 

The analysis of the N = 00 case, carried out in the article, gives an explicit 
analytical illustration of the condition which selects regular solutions. This enables 
to get an insight into the nature of FP solutions of the ERG equations in the LPA. 
Usually these equations are stiff, and gaining a better understanding of the behavior 
of its solutions is quite valuable. 

We would like to mention that the pattern of the curves resembles the pattern of 
eigenvalues (r/ n ) of an eigenvalue problem with a parameter (7). As the parameter 
varies the levels change and various phenomena, like level crossing, may occur. 

The results can be extended to higher order approximations of the Polchinski ERG 
equation. In particular, as it was mentioned in the Introduction, the regular solutions 
studied here can be used as initial input in the iteration procedure used for solving 
the next-to-leading order approximation. Methods of the analysis and the obtained 
results can also be useful in studies of ERG equations of other types in scalar models, 
as well as in other classes of theories, in particular in fermionic theories. 
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